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CLASSIFICATION AND STABILITY OF RELATIVE EQUILIBRIA FOR THE TWO-BODY 
PROBLEM IN THE HYPERBOLIC SPACE OF DIMENSION 2 


LUIS C. GARCIA-NARANJO, JUAN C. MARRERO, ERNESTO PEREZ-CHAVELA, AND MIGUEL RODRfGUEZ-OLMOS 


Abstract. We classify and analyze the stability of all relative equilibria for the two-body problem in the hyperbolic 
space of dimension 2 and we formulate our results in terms of the intrinsic Riemannian data of the problem. 


1. Introduction 


This paper analyzes the relative equilibria of the generalization of the two-body problem to the complete, simply 
connected, hyperbolic, two-dimensional space This is a continuation of the research program initiated in Diacu, 
Perez-Chavela, Santoprete [12], on the study of relative equilibria for the n-body problem in spaces of constant 
curvature. Here we focus on the simplest case of two particles in . We present a rigorous geometric treatment that 
allows us to fully classify and determine the stability properties of all relative equilibria. 


This problem has a long history beginning with the formulation of the Kepler problem in spaces of constant 
curvature. Apparently (see [IHI), the analytic expression of the generalization of Newton’s potential for the Kepler 
problem to the three-dimensional hyperbolic space J-C^ was first given by Schering |3^, following the geometric ideas 
suggested decades earlier in the works of Bolyai [2| and Lobachevsky . The dynamics of this problem is integrable 
and was thoroughly considered by Liebmann |2T|, and more recently by Kozlov et al [T8| and Carihena et al |6]. Just 
like in the euclidean case, the bounded orbits are (hyperbolic) conics and a version of Kepler’s laws holds. A nice 
explanation for this analogy in terms of projective geometrical arguments can be found in [T|. We also mention that 
the ideas behind the generalization of the Kepler problem to constant curvature spaces have been taken further into 
the subriemannian realm m 

The formulation of the n-body problem in spaces of constant curvature follows by generalizing the Kepler potential 
to include the pairwise interaction between all the masses [TB] . In Section 3 of reference [TB] , the authors indicate the 
interest in the investigation of its particular solutions. An explicit expression of the equations of motion can be found 
in |12j . Physical and mathematical motivations to study this problem, as well as some historical details, can be found 
in [TU]. 


We remark that apart from the Kepler and the n-body problem, other mechanical systems with configuration 
manifold a space of constant curvature of dimension two have been discussed very recently, using differential geometric 
tools. This is the case, for instance, of the two centre problem, the harmonic oscillator, and the n vortex problem. 
See [31 m [BJ Uni ISZ] and the references therein. 

Contrary to the euclidean case, the two-body problem in does not reduce to the Kepler problem and is non- 
integrable [ 33 ] • A number of publications have addressed the dynamics of the restricted two-body problem in (see 
null and the references therein). Just like the full two-body problem, it is nonintegrable |23j . 

The impossibility to reduce the two-body problem in Jf^ to the Kepler problem has been associated to the absence 
of “the integral of the center of mass” in [MIIII. From the geometric mechanics perspective, this seems to be related 
to the algebraic differences between the group of orientation preserving isometries of euclidean two-space and of 
(see Remark 4.3). For a treatment of the symmetry reduction of the problem see |5l|33] (see also [35| and the references 
therein). 

Previous results on the existence of relative equilibria for the two-body problem in first appeared in (T^j in the 
case of equal masses and in [14] in the general case. The results in [12] are obtained in the Weierstrass model for 
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J-C^ that arises as an embedding of a hyperboloid in Minkowski’s 3-dimensional space. The authors work with global 
coordinates in the ambient Minkowski space and thus refer to them as extrinsic coordinates. On the other hand, the 
treatment in m is performed in the Poincare disk and upper half-plane models. The authors use the terminology 
intrinsic coordinates and they say that they take an intrinsic approach since these spaces are not embedded in a 
larger ambient space. In this paper we will reserve the use of the term intrinsic to concepts that are intrinsic in the 
Riemannian geometry sense, namely, that are invariant under isometries. In particular, contrary to the conventions 
taken in |I4j . in the terminology of the present paper, a choice of coordinates in a model of or a mathematical 
formula expressed in these coordinates, are not intrinsic. 

Our contribution is to classify all relative equilibria of the problem and to fully describe their stability. Moreover, 
we formulate all of our results in an intrinsic form. In this way, all of our results and formulas (given in section]^ are 
valid in any model of 


We show that the only relative equilibria of the problem arise as a conjugation of the so-called elliptic and hyperbolic 
relative equilibria found in |I4j . In order to study the problem of existence for general configurations, we introduce a 
notion of hyperbolic center of mass between any two masses in Our definition seems to be the correct one from 
a dynamical point of view since it allows us to simplify the study of existence and stability of relative equilibria by 
assuming that the masses lie in a particularly simple configuration. Interestingly, a different definition of hyperbolic 
centroid of mass had been given before in the literature |16j . 

We remark that in the recent paper |30| the authors consider the existence problem of relative equilibria in the 
n-body problem (in particular, in the case n = 2) in the upper half plane of constant negative curvature. In this 
model, they also prove that, up to conjugation, the relative equilibria are elliptic and hyperbolic. However, they do 
not consider general configurations so their description is not exhaustive. In fact, in our paper, we prove that the 
study of existence may be reduced to the case of two masses in canonical configuration (see Proposition 5.4). This 
simplihcation is possible by the introduction of the concept of the hyperbolic center of mass and the exploitation of 
the homogeneity and isotropy of 


Concerning the stability, we show that all hyperbolic relative equilibria are unstable. A similar result for a hyperbolic 
relative equilibrium of three identical masses is given in m- However, we show that the stability of elliptic relative 
equilibria depends on the distance between the masses. If the masses are sufficiently close, the relative equilibrium 
is stable. Otherwise the relative equilibrium is unstable. We give an explicit analytic expression that shows how the 
threshold distance depends on the ratio between the masses and we also show that this distance is related to a critical 
value of the momentum. 


Our stability analysis is based on the Reduced-Energy-Momentum method of Simo et al [36], which gives sufficient 
conditions for a relative equilibrium to be G^^-nonlinearly stable in the sense of Elj. Several comments on this notion 
of stability are in order: First, this is a weaker notion than that of orbital stability, since we regard as G^-stable 
relative equilibria for which nearby initial conditions can drift far away along G^-orbits, the group G^ being the 
subset of the symmetry group fixing the total momentum of the relative equilibrium. It follows from the geometric 
and dynamic nature of Hamiltonian systems with symmetry that to require orbital stability for a Hamiltonian relative 
equilibrium would be too restrictive and as a consequence, G^-stability is the common accepted notion of nonlinear 
stability in this field. An explicit example where a numerical simulation shows a relative equilibrium which is G^- 
stable but not orbitally stable due to this drift can be found in |20j . Second, if a Hamiltonian relative equilibrium is 
G^-stable then the projected equilibrium point to the symplectic reduced space corresponding to its momentum value 
is Lyapounov stable in the usual sense. However the converse is not necessary true since lifting the conditions for 
Lyapounov stability on the reduced space to the original phase space would only take into account initial conditions 
with prescribed momentum. The precise relationship between these two notions of stability is clarihed in [31j . It is 
worth noting that G^-stability of a relative equilibrium also implies Lyapounov stability of the corresponding projected 
point in the reduced Poisson space. However, in this case the relationship between these two concepts of stability is 
more involved. In our case, two conclusions can be drawn about the solutions whose initial conditions are sufficiently 
close to a G^-stable elliptic relative equilibrium. Firstly, they are bounded, and secondly, the variations in the distance 
between the particles is arbitrarily small during the motion. 

Our paper is organized as follows. In section we review known results on the hyperbolic space and of the 
Poincare upper half-plane model H^. In section]^ we formulate the 2-body problem in and we give a useful explicit 
formula for the potential in the H^-model. Section [^describes the symmetries of the system, its infinitesimal versions. 
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and the corresponding conservation laws. In section we introduce the concept of hyperbolic center of mass and we 
prove all of our existence results. Section considers the stability of the relative equilibria found in section Many 
of the results given in sections !-§ are formulated in the H^-model and section is concerned with providing their 
intrinsic form and, as a consequence, the independence of these results on the chosen model. Finally, we include an 
Appendix which contains some standard results on G^-stable relative equilibria for general symmetric Hamiltonian 
systems and for the 2-body problem in an arbitrary Riemannian manifold. 


2. Preliminaries 


By the hyperbolic two-dimensional space we understand the two-dimensional complete simply connected Riemannian 
manifold with constant curvature —1. This is a well defined abstract space that we denote by A well known 
model for it is the upper half plane = {(x, i/) S / y > 0} equipped with the Riemannian metric 


, 2 dx^ + dy^ 
ds = -^- 

yZ 


( 2 . 1 ) 


In this work, we will mainly work with this model but we will formulate our results in such way that they are valid 
in the abstract hyperbolic two-dimensional space Our choice to work in stems from three nice properties. 

The first one is that it can be conveniently covered with a single coordinate chart {x,y). The second one is that this 
model is conformal, in the sense that the angle between any two curves in the upper half plane that intersect is the 
same when measured with the metric and with the euclidean metric. This property suggests the identification 
T(a, that we will use in the sequel. 

The third convenient property is that can be endowed with a Lie group structure arising from the following 
matrix representation: 

()■ 

This Lie group can be interpreted as the semi-direct product of (R^,-) and (K,+). That is, the semi-direct product 
of expansions and translations on M. The Lie algebra I) is spanned by the matrices 


with commutation relation 




0 

0 


[ei, 62] — —Cl- 


A crucial observation is that the hyperbolic metric (2.1) is invariant under the left multiplication on Therefore 
is a subgroup of the group of former isometries of HA This implies that the isometry group of H^ acts transitively 
which means that the hyperbolic two-dimensional space is homogenous. 


The following, well-known result, follows from Hadamard theorem (see, for instance. Chapter 10, Theorem 22 in 

m)- 


Proposition 2.1. 


Given any two points in 'K^ there exists a unique geodesic passing through them. 


The following proposition is also a well-known result (see, for instance. Chapter 3, Example 3.10 in m) 

Proposition 2.2. The geodesic curves on H^ are either half circumferences cutting at right angles the x-axis y = 0, 
or vertical lines x = const. 


Consider a parametrized geodesic ^ : R —>• H^ and let vi,V 2 be tangent vectors to H^ at ^(ti), ^(^ 2 ), for some 
ti,t 2 G R. Suppose moreover, that Vi is perpendicular to £(< 1 ), i = 1,2. We say that the vectors vi and V 2 are equally 
oriented if the bases {vi,£{ti)} and {v 2 , £{£ 2 )} define the same orientation of R^. A pair of vectors that is not equally 
oriented is said to be oppositely oriented. 

Therefore, if £ C H^ is a geodesic in H^, one may introduce the notions of equally and oppositely oriented non-zero 
tangent vectors vi and V 2 based at different points in £ and that are perpendicular to £ (see Figures and for an 
illustration). The correctness of the definition relies on the possibility to parametrize geodesics, and on a choice of 
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(a) Equally oriented. (b) Oppositely oriented. 

Figure 1. Illustration of equally and oppositely oriented vectors for a half circle geodesic in 
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(a) Equally oriented. (b) Oppositely oriented. 

Figure 2. Illustration of equally and oppositely oriented vectors for a vertical geodesic in 


orientation of Hence, this concept may be transferred, in a natural way, to the abstract hyperbolic space 
This notion will be useful in the intrinsic description of relative equilibria in section [Tj 

For future reference, we note that the Riemannian distance, d, between two points {xi,yi), {x2,y2) G satisfies 

C0.h(d) = 1 + + (K - y^)\ (2.2) 

2yi?/2 

In particular, if {xi,yi) = (±cos0,sin0), with 0 < 9 < |, and {x 2 ,y 2 ) = (0,1) then 

sind = sech(d), cos 0 = tanh((i). (2.3) 


2.1. The group of isometries of Let 

SL(2, K) = {5 e GL(2, K) | det 5 = 1}, 

where GL(2,K) is the group of invertible 2x2 real matrices, be the special linear real group, which is a 3-dimensional 
simply connected, real Lie group. 

The group SL(2,K) defines an action on as follows. Given a matrix g G SL(2,M), 



and an element (x,y) G H^, we define ^'(g, (x,y)) via the Mobius transformation 

a{x + iy) +b 
c{x + iy) + d’ 


'i>{g,{x + iy)) 


(2.4) 
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where we have used complex notation for elements in One can check that the upper half plane y > 0 is invariant 
under the above transformation and that the action axioms are verified, so is well defined. 

It is well known that for any g G SL(2,K), the map '^{g, •) is an orientation preserving Riemannian isometry of 
equipped with the hyperbolic metric. Moreover, the group of orientation preserving, proper isometries of is the 

quotient group SL(2,]R)/{±/} PSL(2,K). The quotient is taken to account for the fact that '^{g, •) = •)• Iii 

the sequel we will sometimes work with the group SL(2,K) keeping in mind that isometries are counted twice. 

As mentioned before, a good number of isometries defined by elements in SL(2, K) correspond to left translation on 
H^. Indeed, left multiplication by the element 


f yo xo 

I 0 1 


e 


corresponds to the action of the matrix 

i I e SL(2,K). 

^ VvH J 

In fact, we can consider as a Lie subgroup of SL(2,M) via the group homomorphism 


{ yo a;o 
1 0 I 






(2.5) 


The action ik of SL(2,IR) on introduces a new type of symmetries that do not arise as translations on 
These transformations are analogous to usual rotations on the plane. This family of symmetries is generated by an 
subgroup of PSL(2,IR) consisting of matrices of the form 


9t ■■= 



— sin 
cos ( 


) 


t S M mod 2?!. 


The corresponding transformation on acts like a rotation by the angle t around (0,1) in the following way: it fixes 
the identity element (cc = 0,?/ = I) and takes the geodesic through (0,1) with tangent vector i; € = r(o,i)IHl^ into 

the the geodesic through (0,1) with tangent vector g- 2 tv G = T(o.i)IHI^ . Thus, this subgroup acts transitively 
in the space of directions through (0,1). This, together with the homogeneity of space, shows that the hyperbolic 
two-dimensional space is isotropic. 

We shall call the above transformations “hyperbolic rotations” around (0,1). 

It is easy to prove that the hyperbolic rotations around (0,1) are just the isotropy subgroup SL(2,K)jp of (0,1) 
with respect to the SL(2, M)-action on Thus, using Theorem A.I.l in [19], the compactness of SL(2 ,IR)(q and 
the fact that the action of SL(2,IR) on is transitive, we deduce an important result: the action of SL(2,IR) on 
is proper. 


2.2. A property of isometries in In this section, we present a basic result about isometries of that will be 
useful in the sequel. 

Proposition 2.3. Let i C be a complete geodesic curve and let {xQ,yo) G i. There exists an isometry o/H^ that 
maps I onto the upper half of the unit circle and (xo,yo) onto (0,1). 

Proof. By homogeneity of there exists an isometry defined by gi G SL(2,IR) that maps (xo,yo) onto (0,1). Since 
isometries take geodesics into geodesics, the image of £ by this isometry is a certain geodesic £ through (0,1). Next, 
by isotropy of there exists a hyperbolic rotation around (0,1) defined by 52 G SL(2, K) that maps £ onto the upper 
half of the unit circle. The element 5251 € SL(2,K) defines the desired isometry. □ 

An explicit formula can be given for the isometry described in the statement of the above proposition. Let v G 
T(^XQ,yo)^ be any non-zero vector tangent to £. It can be written as 

V = A2/o(cos6*o,sin0o) G 
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where 


1 




COS I 


^ cos(^)a:o+s 


1 . f0o\ 

sin ( y ) , 


vyo \ - / 

satisfies the requirements of Proposition |2.3| 

Note that go = < 7251 , with gi,g 2 G SL(2,IR) given by 

/ 1 _J£o_ \ 

51 = ^yg 2 


- sin (^) CCO + cos (|i) j /0 

7^^ ’ 


cos(|i: 

sin(^) 


-sinm 

cos(^) 


3. Geometric formulation of the 2-body problem in 

The 2-body problem in can be defined as the simple mechanical system on the configuration space Q = A, 

where A is the set of collision configurations, and whose Lagrangian L : TQ —M is given by 

L = i {mivl + TO 2 v|) - V{d). 

In the above equation mi and vf are, respectively, the mass and the square of the hyperbolic norm of the velocity 
of the particle. The positive number d is the hyperbolic distance between the first and the second mass and the 
potential V : M"*" —M is given by 

V(d) = —kmim 2 coth(d), (3.1) 

where fc is a positive constant. The hyperbolic cotangent dependence of V on the distance is well-accepted to be the 
“correct” generalization of the Newtonian gravitational potential for two reasons. The first one is that the Kepler 
problem with this potential satisfies Bertrand’s property, namely, all of the bounded orbits are closed. The second 
reason is that, just like the Newtonian potential Vff{d) = ^ is (proportional to) the fundamental solution of the 
Laplacian operator on K^, the proposed potential V{d) is (proportional to) the fundamental solution of the Laplace- 
Beltrami operator on As a consequence, it satisfies the hyperbolic version of Gauss law. 

In what follows we will work in the model. We will also denote Q = (]3^)^\ (A) with elements q = {xi,yi,X 2 , j/ 2 )- 
The Hamiltonian of the problem is h : T* Q ^ M. given by 

HPq) = \ \\Pq\\^ + y{<l) 

where || • || denotes the norm in the fibers of T*Q induced by the kinetic energy Riemannian metric on Q given by 


dx^ -I- dyi 

mi - 2 - 

vi 


7712- 


dx\ 


dyl 


Vl 


(3.2) 


Using (2.2) 




2 yi ?/2 


and the relationship 


coth(cosh ^(z)) = 


that is valid for all z > 1 we can write 


'J — 1 

{xi - X2f + yl + yl 


7 {{xi - X 2 y + (yi - y 2 )^)((a;i - 3^2)^ + (yi + 2/2)^) 


V{xi,yi,X2,y2) = -kmim2 


(3.3) 
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4. Infinitesimal symmetries and the momentum map 


Since the potential energy only depends on the hyperbolic distance between the particles, it follows that the diagonal 
action of SL(2,IR) on Q lifts to a symmetry of the problem. It is well known that nontrivial isometries of that fix 
more than one point do not preserve orientation. Therefore, since we have removed the collision configurations, this 
action is free. 

In addition, using the proper character of the action of SL(2,IR) on H^, we directly deduce that the action of 
SL(2,IR) on Q also is proper. 


Remark 4.1. Since the action of SL(2,K) on Q is free and proper, we have that Q is the total space of a principal 
SL(2, K)-bundle over the space of orbits (5/SL(2,R). In fact, using Proposition 2.3 and formula (2.3), one may prove 
that (5/SL(2,IR) is diffeomorphic to K+ = {r S M : r > 0}. Moreover, under this diffeomorphism, the orbit projection 
TT : <5 M’*' is just the restriction to Q of the hyperbolic Riemannian distance. o 


In this section we write the inhnitesimal version of the SL(2, R)-symmetry in the model, we make a classification 
of the Lie algebra elements that will be useful for the sequel and we compute the corresponding momentum map. 
The Lie algebra of both PSL(2,R) and SL(2,IR) is the 3-dimensional real linear space 

s[(2,K) = {^ e M2x 2(]R) I traced = 0 }. 


Hence, any non-zero element in sl(2, K) is of one of the following three types: 

(i) Elliptic. These elements have two complex conjugate, purely imaginary eigenvalues. 

(ii) Hyperbolic. They posses two real eigenvalues with the same absolute value and 
opposite signs. 

(iii) Parabolic. They have a multiplicity two zero eigenvalue (and are not diagonalizable). 


The following are representatives of elements of the above types and form a basis of s[(2,R). 


^e:= 


? 

4 0 


■ = 


0 


0 -T 


ip 


We have that 

[^ei C/i] — Ce + Cpj Cp- (4-1) 

Let i be an arbitrary elliptic element in sI(2,IR) with eigenvalues w G R, and corresponding complex eigen¬ 

vectors uztiv with M, n € R^. We then have 


iu = --V, 


c 

iv = -u. 


(4.2) 


From these relations it is easy to show that u and v are linearly independent over R^. Suppose that {u, n} is positively 
oriented. Then one can scale u and v by the same positive factor so that ( |4.2[ ) continues to hold and the matrix 
g =: (u|'c), that has u and v as columns, has determinant one so it belongs to SL(2,R). In view of (4.2) we obtain 

Adg-i^ = g~'^ig = -w^e- 

If {u, n} is not positively oriented, after a rescaling by a positive factor, the matrix g = {u\ — v) £ SL(2, R) and satisfies 
Adg-i^ = In a similar way, given a parabolic element i G s[(2,R), one can find a basis of R^ such that 

iu = 0, iv = u. 


If {u, n} is positively oriented then we can scale u and v by the same positive factor so that the determinant of the 
matrix g = (u|n) is I, and the above relations still hold. In this case 

= ip. 

If {u,v} is not positively oriented one can rescale by a positive factor so that the matrix g = {u\ — v) G SL(2,R) and 
we have Adg-i^ = —ip. 

Finally, any hyperbolic element in ^ G s[(2,R) with non-zero real eigenvalues is diagonalizable. Rescaling the 
eigenvectors (by positive or negative factors) shows that one may choose g G SL(2,R) such that 

i = ^9ihg~^ = i^Mgie. 
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In conclusion, we have proved the following. 


Proposition 4.2. Let ^ be an element ms[(2,K). Then: 

(i) ^ is elliptic (respectively, hyperbolie) if and only if there exists g € SL(2,K) and w € K, w ^ 0 such that 

5 = wAdg^e 

(respectively, ^ = ujKdgf^h)- 

(ii) ^ is parabolic if and only if there exists g S SL(2,IR) such that 

C = ±Adg^p. 


The description of s((2,IR) in terms of elliptic, hyperbolic and parabolic elements ties in very well with our study 
of SL(2,IR) done in section |2.1| By a direct calculation, one can see that elements and fp span a two-dimensional 


(4.3) 


Lie subalgebra of s((2,K). The underlying Lie subgroup of SL(2,]R) is isomorphic to via the map (2.5). 
Starting from the definition of the action in (2.4), a direct calculation shows that 

^ (exp(a;^/,t), {x, y)) = e“*(a:, y), T (exp(wCpt), {x, y)) = {x + cot, y). 

On the other hand, one has 

(sin(a;t)(a;^ -I- — 1) -I- 2a:cos(a;t), 2y) 


4' (exp(a;^et), {x,y)) = 


+ y'^ + 1) -\- {1 — x'^ — y'^) cos(a;t) -I- 2x sin(wt) ’ 


(4.4) 


which is a hyperbolic rotation around (x,y) as introduced in Section 2.1 The terminology is somewhat unfortunate 
since “elliptic” elements in the Lie algebra sl(2, K) generate “hyperbolic” rotations. 

Putting a; = 1, differentiating and evaluating at t = 0 formulas (4.3) and (4.4), we obtain the following expressions 
for the inhnitesimal generators of the action 

{£.h)m^ix,y) = {x,y) G i^p)M^{x,y) = (1,0) G , 

— a;^ — 1 


i^e)n^{x,y) = -1- -,-xy'^ G T(^x,y)^ 


The corresponding expressions for the infinitesimal generators of the diagonal action of SL(2,I 
obtained from the above formulas. 

In our treatment we will identify the dual Lie algebra s[(2,]R)* with s[(2,]R) via the pairing 
(Ai,0 = 2Trace(/rC), for y G sl(2,M)* ^ s[(2, K), ^ G sl(2,M). 

The basis of s[(2,IR)* = s[(2,M) that is dual to via this pairing is {ph, Te, Tp} with 

h-h — Me = Cp) Mp “ (?p T Ce)- 

In addition, with the previous identification 

ad|/r = [/r,^], and Adg-iy = yyy“^ 


(4.5) 

on Q are readily 

(4.6) 

(4.7) 

(4.8) 


We can now compute the momentum map associated with this symmetry. The momentum map of the cotangent 
lift of the diagonal action of SL(2,IR) on Q is J : T*Q —>■ s[(2,IR)* = s[(2,IR) given by 

^ Pxi {yf - xl - 1) - 2py.x^yi ' 


Jixi,X2,yi,y2,Pxi,Px2,Pyi,Py:,) = '^XiPxi + ViPy^ + X! 


Me 


(4.9) 


+ ] pxij Mp- 


The components of the above expression can be checked to be integrals of the equations of motion. Since the 
symmetries associated to and are related to the homogeneity of space, the components of ph and Pp are natural 
generalizations of linear momentum. On the other hand, since the symmetry defined by is a consequence of isotropy 
of the space, the component of pe is naturally analogous to the classical angular momentum. Our interpretation of 
the integrals of motion differs from that of m- 
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Remark 4.3. The system has the same number of conserved quantities as the euclidean two-body problem arising 
from the action of the group of orientation preserving isometries. However notice that the subgroup representing 
translations is in the euclidean case whereas for the hyperbolic case it is A fundamental difference between 
the two is that the latter is not abelian. As a consequence, the symplectic reduction of T*Q by yields a reduced 
space whose generic dimension is 6 (the stabilizer of the momentum is generically trivial). This contrasts with the 
euclidean case where the symplectic reduction by the action of has dimension 4 and the resulting reduced system is 
equivalent to the euclidean Kepler problem in the plane. The two extra dimensions that are lowered in the reduction 
of the euclidean case are interpreted as passing to the center of mass coordinates, a procedure that is not possible in 
the hyperbolic case. 

Another difference is related to the fact that is a normal subgroup of the euclidean group SE(2,M) so one can 
perform reduction by stages [5S] and further reduce the system by rotations. The ultimately reduced system is the 
Hamiltonian system with one-degree of freedom for the radius under the influence of the effective potential of the 
Kepler problem. On the other hand, is not a normal subgroup of SL(2,K) and, as a consequence, the hypotheses 
for reduction by stages [5S] are not met. 

o 


5. Existence of relative equilibria 

In this section we will classify all relative equilibria of the problem. We will show that the only relative equilibria 
arise as conjugation of the ones found in [14]. 

We begin with the following. 

Definition 5.1. The hyperbolic center of mass of two masses mi, m 2 in is the unique point in that satisfies 
the following properties: 

(i) It lies along the (unique) geodesic connecting mi and m 2 . 

(ii) If di denotes its distance to mi and d 2 its distance to m 2 , then 

mi sinh(2di) = m 2 sinh( 2 (i 2 ). 

Note that in the above definition di + d 2 is the distance between mi and m 2 . It can be checked that the hyperbolic 
center of mass of two masses is well and uniquely defined. Figure [^illustrates our definition. 


m 


1 


Figure 3. Schematic representation of the hyperbolic center of mass between particles mi and m 2 . The 
curve connecting the particles is a geodesic and the relation mi sinh(2di) = m 2 sinh( 2 d 2 ) is satisfied. 



Suppose that Xi,X 2 G denote the positions of the masses mi and m 2 and that C G is their hyperbolic 
center of mass. It is clear that the hyperbolic center of mass of the masses mi and m 2 at the positions g ■ Xi,g ■ X 2 
is g ■ C where g denotes any isometry of 
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Our definition of the Iwperbolic center of mass is the one appropriate for our purposes of classifying the relative 
equilibria of the problemjj Interestingly, our definition differs from the position of the hyperbolic centroid defined in 
m by generalizing the euclidean “lever rule”. In this treatment, the hyperbolic centroid is a material point with a 
certain mass. For us, the hyperbolic center of mass is just an abstract point in 5f^. 

Definition 5.2. Two masses mi, m 2 in are said to be in canonical configuration if they are respectively located at 
(cos 6 * 1 ,sin 6 * 1 ) and (—cos 02 jsin 6 * 2 ), where the angles 0 i ,02 satisfy 0 < 0 i ,02 < f > o-nd the following relation holds 

mi cos02sin^0i , _ 

— = —9 -. (5.1) 

m 2 sin 62 cos 01 

The definition for a canonical configuration is illustrated in Figure]^ 



Figure 4. Two particles of masses mi and m 2 in canonical configuration. The relation (5.1 1 holds. 


Using (2.3), we deduce the following characterization of two masses in canonical configuration. 


Lemma 5.3. Two masses in which are located at (cos0i, sin0i) and (—cos 02 , sin 02 ), with 0 < 0 i ,02 < f, are in 
canonical configuration if and only if its center of mass is the point ( 0 , 1 ). 

Now, we may prove the following result 

Proposition 5.4. Given two masses mi, m 2 in there exists an element o/SL(2, K) that maps them to the canonical 
configuration. 


Proof. In view of Proposition 2.3 there exists an isometry that maps the geodesic passing through mi and m 2 onto the 
upper unit circle with the center of mass mapping onto (0,1). If after this transformation mi lies in the first quadrant 
we are done. Otherwise we apply the isometry defined by the matrix 

0 -1 


1 0 


e SL(2,] 


□ 


Relative equilibria are solutions of the equations of motion that are contained in a group orbit for the action of 
SL(2, K) on T*Q. It is well known that a relative equilibrium is equivalent to a pair {q, € Q x q, and in this case the 


solution is given by z{t) = exp(t^) • {{f.Qiq),-)), where ((•, •)) denotes the Riemannian metric (3.2) on Q. The element 
^ G g is called the velocity of the relative equilibrium. We will refer interchangeably to both z{t) and Pq = {{£,Q{q), •)) 
as the relative equilibrium. 

If z(t) is a relative equilibrium then so is 5 • z(t) for any g € G with the same stability properties as z(t). The 
velocity of g ■ z{t) is Adg^. 


^We are not claiming any general properties of the evolution of the hyperbolic center of mass under the dynamics, see Remark 7.3 
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Due to the equivariance of the momentum map, a relative equilibrium pq with velocity ^ and J{pq) = p must satisfy 

ad^/i = 0. (5.2) 

Pairs (q, corresponding to relative equilibria are characterized by the condition 

dV^{q)=0, (5.3) 

where Vj S C°°{Q) is the augmented potential defined by 


Here I{q) : sl{2, 


V^iq) = V{q)--{I{q)^,0. 

sl(2,IR)* is the locked inertia tensor defined as 


(5.4) 


Since relative equilibria come in group orbits, according to Proposition |5.4[ we restrict our study of existence of 
relative equilibria assuming that the masses are in canonical configuration. 


Let Pq be a relative equilibrium with velocity ^ + P^p- One can readily check, using (4.91 that 

H 

P — cos 9 1 cos 02 E 


( m2(cos02 + COS0i) / iJ (1 — 2cos0i COS02)P +COS01 COS02L^ 

d\Pq)— I P _ r-no fl, P —}J 


2 sin 02 cos 01 


(5.5) 


where we have used (5.1). 


We start the classification of relative equilibria by applying condition ( |5.2[ ) in order to isolate possible candidates. 
In view of our identification (4.6) this is equivalent to the condition [ J {pq), ^] = 0 where [•, •] is the matrix commutator. 


Up to an unessential non-vanishing factor this commutator is 

P{E — P) —E[{{1 + cos 01 cos02)i? — (1 + 2 cos 01 COS02)P) 

-iL((l-kcos0iCOS02)L;-P) -P{E-P) 

There are only three non-trivial possibilities for the above matrix to be equal to zero: 

(a) E = P = 0, E[ = uj^0. The Lie algebra generator ^ = uj^h is hyperbolic. 

(b) H = P = 0, E = uj^0. The Lie algebra generator ^ is elliptic. 

(c) H = 0, P = E = uj^0, The Lie algebra generator ^ = a;(^e + ^p) is hyperbolic. 

In particular notice that the above conditions do not include the possibility of having H = E = 0 and P ^ 0. This, 


using Proposition 4.2 and 5.4 gives an alternative proof to the following result given before in El- 


Proposition 5.5. There do not exist relative equilibria for the two-body problem in with a parabolic Lie algebra 
generator. 


We now use equation (5.3) to investigate when possibilities (a), (b), (c), are realized for the particular potential 
energy ([O. 


and {/Tg,/x/i,/Ip}. 


Using (|4.5|) and (|5.4|) we obtain the following representation for the locked inertia tensor referred to the bases 

I{xi,yi,X 2 ,y 2 ) = 




1 miiixf + yf + 1)2 - 4y2) 


Z=1 


yt 


/ll 

h2 

Il 3 \ 


h2 

I22 

I23 , 


Il 3 

I23 

I33 J 


2 

y^mi 

2 = 1 

ix^ + yf) 

yf ’ 

2 

^33 = 2^ ^ 

i=i 


1 xfyf) 


2=1 


Vi 


T _ 1 + - 2 / 1 ) 

-'13 -~o2^ - -b2- 


hs - 


rUiXi 


2 = 1 




with 

















12 


L. C. GARCfA-NARANJO, J. C. MARRERO, E. PEREZ-CHAVELA, AND M. RODRIGUEZ-OLMOS 


With the above formulae in hand, one obtains the following expressions for the augmented potential Vj in cases 
(a), (b) and (c). 


Vui^{xi,yx,X2,y2) = V{xi,yx,X2,y2) 
Vu,^^{xi,yi,X2,y2) = V{xi,yi,X2,y2) 
^^^(«e+«p)( 2 ;i-yi> 2 ^ 2 , 2 / 2 ) = V{xi,yi,X 2 ,y 2 ) 



2=1 


mtix'l + yi) 

yf ’ 

m^{x^ + ivi - iy){x^ + iVi + 1)^) 

2 ’ 
Vi 

mi((xi - ly + yiYyixi + ly + y1) 
0,2 


A direct calculation with MAPLE'^'^ that is too long to be included here, shows that, in cases (a) and (b), 
q = (cos 01 , sin 01 , —cos 02 , sin 6 * 2 ) is a critical point of if and only if 

2 sin^ 02 sin^ 0i / sin^ ^2 \ . sin^ 02 sin^ 0i / sin^ 0i \ 

^ (cos 02 + cos 01 \COS 02 ”^^/ (cos 02 + cos 01)2 \ COS 0i 

On the other hand, if ^ satisfies condition (c), then a calculation with MAPLE"'"^ shows that condition 

mlk cos 02 sin^ 62 sin^ 0 i (cos 0 i + cos 62 ) = 0 . 

Considering that k,m 2 > 0, and the restrictions on 0 i ,02 coming from Definition |5.2[ it follows that this condition 
can never hold. Therefore we have shown the following. 


5.3) implies 


Theorem 5.6. For any values of mi, m 2 > 0, there exist exactly two families of relative equilibria for the two body 
problem in H^. Representatives of these families in a canonical configuration q G Q are {q,ujfh) and {q,ujfe) where lu 
satisfies (5.6). 


The existence of a representative of the relative equilibrium {q,u}fh) in the general case mi y rn 2 is proved in 
Theorem 7 of [14]. The condition of this theorem, given in equation (96|^ is readily seen to be equivalent to (5.1). 
The existence of the relative equilibrium {q,u}fe) for mi y m 2 was established in Theorem 3 of [14j using the Poincare 
disk model. We have just shown that, up to conjugation, these are the only ones. Following the terminology in 
and in accordance to the nature of the Lie algebra generator, we shall refer to these families as hyperbolic or elliptic, 
respectively. 

Notice that the orbits in generated by hyperbolic (respectively elliptic) generators are unbounded (respectively 
bounded). The same is true for hyperbolic and elliptic relative equilibria in Q. 


Remark 5.7. In chapter 4 of |9] a relative equilibrium is called elliptic if the one-parameter subgroup {exp(,^t) : f S K} 
generated by its velocity f, is dense in a torus subgroup of the corresponding symmetry group. Relative equilibria that 
are not elliptic are called runaway relative equilibria. In our case, the uni-parametric subgroup {exp(^eO : t S K} of 
SL(2, K) is closed and isomorphic to (a one-dimensional torus) while the uni-parametric group generated by fh is not 
compact (and hence its closure cannot be contained in a torus). Thus, our terminology for elliptic relative equilibria 
is consistent with that of |9], while our hyperbolic relative equilibria correspond to runaway relative equilibria. From 
the general theory developed in |5] it follows that our elliptic relative equilibria are periodic whereas our hyperbolic 
relative equilibria “run out” of any compact set K C T*Q as t —>■ ±00 (i.e. they are unbounded curves). o 


6. Nonlinear stability 

Let z : K —> T*Q be one of the relative equilibria found in Section and J(z(0)) = p, the initial momentum 
associated with z. Then, as we know, J{z{t)) = p, for all t. So, if z{t) = exp(^t) • z(0), with ^ € 0 then, using 
the equivariance of the momentum map J, we deduce that p = J{exp{ft) ■ z(0)) = exp(^t) • p, which implies that 
exp(,^t) belongs to G^, where is the stabilizer of the momentum p under the coadjoint representation of SL(2,M) 
on s[(2,IR)*. Thus, the relative equilibrium z is contained in the G^-orbit through z(0). 


^We note that there is a typo in this equation and it should read as it is easily deduced from the proof of the theorem. 

J sr ^2 P2y2 ^ 
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In this section we discuss the G^-stability of the relative equilibrium z{t). The precise definition of G^-stability 
(see [31]) is the following: 

A relative equilibrium z : M. ^ T*Q is said to be G^-stable if for every G^-invariant open neighborhood U of z(IR) 
there exists an open subset W in T*Q, such that z(M) C IT C [/, and with the property that for every solution z' 
with initial condition in W, we have that z'(t) G U, for all t . 

Since we are dealing with a simple mechanical system, we use the reduced Energy-Momentum method (REM) of 
Simo et al [36] . which studies the signature of the restricted Hessian Here, is the augmented Hamiltonian 

defined as (see | 21 |) 

and iV is a complement to • Pq in T 2 J“^(/i), where p = J(pq)- 

The reduced Energy-Momentum provides a convenient block-diagonalization of this bilinear form. Let Pq be a 
relative equilibrium with velocity ^ and momentum p. Then 

Ar 0 0 

0 {d^V^ + cott)\v,^^ 0 ) ( 6 . 1 ) 

0 0 a: 




relative to a splitting 


N ~ Trig © V)„t © V: 


( 6 . 2 ) 

Here AT is a positive definite bilinear form on Id*^ and equations (6.1) and ( |6.2| ) (and therefore, the reduced Energy- 
Momentum method) hold if the bilinear form Ar : td-ig x V)ig —)■ K is non-degenerate. 

Assuming that G^ is compact, a sufficient condition for G^-stability of the relative equilibrium pq is that d'^^h^\N be 
definite. If the conditions of REM hold, then this is equivalent to both Ar and + corr being positive definite. On 
the other hand, an odd number of negative eigenvalues of dl^h\N implies spectral instability and therefore, instability 
regardless of the compacity of G^ | 21 |. 


Now we define the different spaces appearing in the splitting (6.2). First, 

Tig = {A e g : {I{q)X, 7 ]) = 0 V ?7 e g,,} 

where g^^ is the Lie algebra of Gfj,. Now let (g^^ • q)-^ C TqQ be the orthogonal complement of g^^ • 9 with respect to the 
Riemannian metric (3.2), then 


Mnt = {u e (g,, • q)-^ : I{q) ^{{DI ■ u)(^)) S g^ } . 


(6.3) 


The block Ar is defined by 


Ar(Ai, A 2 ) = (ad^^/r, J(g) + a.<lx^{I{q) V)) , 


while 


corr(ui,U 2 ) = {{DI ■ Vi ){^), I{q) \{DI ■ V2)i0)) ■ 

For the sequel, it is useful to notice that in a canonical configuration, using the expression for the locked 

inertia tensor simplifies to 

cos 62 

0 

— cos 02 


,02) = 


m 2 (cos 01 + cos 02 ) 

sin^ 09 



(6.4) 


6.1. Hyperbolic case. We consider a hyperbolic relative equilibrium (q,^) where q = (cos0i,sin0i,— cos 02 ,sin 02 ), 
^ = uj^h, and where condition (5.6) holds. Using (5.5) we obtain 

(cos 01 


p = J{pq) = ujm2- 


cos 02 ) 
sin^ 02 cos 01 


-Ph- 


where uj satisfies (5.6) 


In view of (4.1), (4.7) and (4.8) we obtain that g^ is spanned by ^h- 
Using (6.4) the space Tig is spanned by 


Cl = Ce, 


C2 = ip- 
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With respect to this basis, the bilinear form Ar is represented by the matrix 


Ar = 


(cos 01 + cos 02 ) cos 02 f 1 
sin^ 02(1 — cos 01 cos 02 ) 


-1 


-1 

_ 1 

COS^ 9i COS^ 02 


Since Arn and det(Ar) are both positive, then Ar is positive definite and therefore REM is applicable since (6.1) 
and (| 6 . 2 |) hold. 


We now compute the internal space Idnt- Since is spanned by ^h, in view of (4.5) we obtain 

Qfi - q = span {(cos 01 , sin 0 i, — cos 02 , sin 02 )} . 

From this 


( 0 M • 9)-^ = <^ - 


sin 01 


b- 


sin 02 


cos 01 cos 02 


a,b,c] : a, b,c G 


If we now use (6.3) we find that the space I^nt is one-dimensional. After choosing a particular generator 
V = [— sin^ 01 cos0i(cos^ 02 - 1 - 1 ), sin0i cos^ 0i(cos^ 02 -I- 1), sin^ 02 cos 02 (cos^ 0i -I- 1), sin 02 cos^ 02 (cos^ 0i -I- 1)) , 
we obtain 


(dgV} + corr)|y.„^(u,u) = - 


km 2 cos 6*2 sin^ ^i(cos^i cos 6*2 + l)(cos^ 6*1 + sin^ 61 cos^ 62 + 3) 
(cos 61 + cos 6 * 2 ) cos(^i) 


(6.5) 


which is negative. 


Therefore, since dp has signature (+, +,—,+), then every hyperbolic relative equilibrium is unstable. 


6.2. Elliptic case. We consider an elliptic relative equilibrium (g, ^) where q = (cos0i, sin0i, — cos 02 , sin 02 ), ^ 
and where condition (5.1) holds. Using ( |4.9| ) we obtain 

. COS02(COS01 + COS02)w , , 

M = JiPq) = m 2 -r^---(/ie - Up). 


In view of (4.7) we obtain that 


fj. = —m2- 


sin^ 02 

cos 02 (cos 01 + COS 02 )w 
sin^ 02 


Ce. 


( 6 . 6 ) 


From this is clear that Gp ~ with Lie algebra = spanj^ej. In particular, Gp is compact. 
We now compute l}ig. From expression (6.4) we find that it is spanned by 

Cl — Chi C 2 = Ce + Cp- 

With respect to this basis, the bilinear form Ar is represented by the matrix 


_ ^ 20 ^^ (cos 01 + cos 02 ) cos 02 / 
sin^ 02 V 


COS 02 


0 

1 + cos 01 cos 02 


Since Arn and Ar 22 are both positive, then Ar is positive definite and therefore REM is applicable since (6.1) and 
( 1 ^ hold. 


We now compute the internal space l}nt- We start by obtaining Qp ■ q. From the above expression for Qp, and (4.5) 
we obtain 

Qp • q = span | (— cos^ 0 i, — cos 0 i sin 0 i, — cos^ 02 , cos 02 sin 02 )} • 

From this 


(flp ■ = 


c sin 02 — 0 cos 01 — 6 cos 02 , , 

a, -—- ,b,c] : a,b,cG 

sin 01 


If we now use (6.3) we find that the space Vint is one-dimensional. After choosing a particular generator 
w = [— sin^ 01 cos0i(cos^ 02 -I- 1), sin0i cos^ 0i(cos^ 02 -I- I), sin^ 02 cos 02 (cos^ 0i + 1), sin 02 cos^ 02 (cos^ 0i + I)) 
we obtain 

m|fcu(I — '«^)^(1 -I- uv)F{u, v) 


(d^U^ +corr)|y„,(w,u;) = 


i(tt -I- v) 


(6.7) 
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with 


F{u, v) = 1 — — U‘‘ — V 


where u = cos0i and v = cos 02 - It is readily seen that the sign of (6.7) is the sign of F(u,v). 
Notice that u and v are related by (5.1). We can express 

— 1 + yj(u^ — 1)2 -I- 4c^u^ 


j(u) = 


2 cu 


( 6 . 8 ) 


(6.9) 


where c = —. It is readily seen that 

m.o 


lim v(u) = 0 , lim v(u) = 1 , v'(u) > 0 . 

It —>-0 li —>-1 

Therefore, (d^V^ + corr)|v:_^j(i(;,i(;) is a positive multiple of f(u) = F(u,v(u)). It is easy to see that 


lim f(u) = 1, 

w —>-0 


lim f(u) = —4. 


Using the chain rule we find that / is decreasing and therefore f(u) changes sign at a unique value uq G (0,1). 
Therefore, for 0 < u < mq the signature of is (+,+,+,+) and the elliptic relative equilibrium is nonlinearly 

stable, and for mq < u < 1 the signature of dp^h^j/^/ is (+, +, —, +) and the elliptic relative equilibrium is unstable. 

Substitution of (6.9) into ( 6 . 8 ) shows, after some algebraic manipulations, that uq is characterized as the unique 
zero between 0 and 1 of the polynomial 

p(x) = 3a;® + (16c^ — 8 )a;® + 6 a;"^ — 1 = (3a;^ + l)(a;^ — 1)® + 16c^a;®. 

More precisely, the elliptic relative equilibrium is stable if p(u) < 0 and unstable if p(u) > 0. Recalling that u = cos0i, 
the condition for stability can therefore be written as 


nil sin® 9iy/3 cos^ 0 i + 1 
TO 2 4 cos^ 01 


( 6 . 10 ) 


This formula will be useful to give an intrinsic description of the stability conditions of elliptic relative equilibria in 
section 0 


In view of (6.9), the momentum ( 6 . 6 ) can be written as a function of u. The change in stability is related to a critical 
value of the momentum fj.(u) as the following proposition shows. This behavior of the stability of a parametrized family 
of relative equilibria with respect to a critical point of the norm of the momentum seems to be more than a mere 
coincidence since it has been observed in a similar context in other systems. See for instance Remark 3.13 in [28) . 

Proposition 6.1. Let || • || be any vector space norm in s[(2,K)*; then ||/i(M)|| is increasing for 0 < u < uq and 
decreasing for uq < u < 1, where uq is the unique zero of f{u) = F{u,v{u)) in the interval (0,1). 


Proof. Substitution of (5.6) into ( 6 . 6 ) (using ( |5.1[ )) yields 

pL = =F\/fcTO2mi\/cOS0i(l — COS^ 02)'fe = p\/km 2 miyju{l — 

Therefore, 

||/i(u)|| = Aa/u(1 - v"^), 

where the positive number A = \/krri 2 rni\\^e\\- We have 

d||^(M)|| l — v'^ — 2uvv' 

= A - 


du 


2 ^u{l — n^) 


On the other hand, the relation (5.1) is written in terms of u and v as 

1 — 


1 — 


Implicit differentiation with respect to u yields 


(u^ + 1)(1 — v"^)v 
(1 — M^)(l + v'^)u' 


V 
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Hence 

^{1 — v^){1 — ^ {1 — v‘^)F{u,v) 

du 2yju{\ — — it2)(l + ti2) 2^Ju{l — i;2)(l — u'^)[l + v"^) ’ 

that has the sign of F{u,v). 


□ 


Two physical implications of G^-stability may be given using a general result for symmetric Hamiltonian systems 
and a standard result for the 2-body problem in an arbitrary Riemannian manifold (see the Appendix). In fact, let 
z : K —>■ T*Q be a G^-stable relative equilibrium. Then, we have: 

(i) All the solutions sufficiently close to z are bounded. 

(ii) If ro > 0 is the constant hyperbolic distance between the two particles along z, it follows that the hyperbolic 
distance between the two particles remains arbitrarily close to rg along all the solutions that are sufficiently 
close to z. 


Properties (i) and (ii) follow using Propositions A.l and A.2 (in the Appendix), respectively. For (i) to hold it is 
essential that G^ is compact. 


7. Intrinsic description 

We now give an intrinsic description of the relative equilibria of the problem involving only Riemannian data and 
without referring to a canonical configuration. Therefore, this description will be valid for any model of 

Theorem 7.1. Consider two masses mi and m 2 in the abstract hyperbolie two-dimensional spaee TC^. A hyperbolie 
relative equilibrium corresponds to the solution characterized in the following way. Let di, i = 1^2 denote the distance 
between the mass mi and the hyperbolie center of mass of mi and m 2 , d := di d 2 the distanee between the masses, 
and let oP satisfy: 

2 2kmi 2km2 

sinh^ ( d) sinh ( 2 ^ 2 ) sinh^ ( d) sinh (2di) 

Then, at every instant during the motion, 

(i) di is constant, 

(ii) the velocity vectors vi and V 2 are perpendicular to the unique geodesic that contains mi and m 2 and are 
equally oriented, 

(iii) the hyperbolie norm of Vi is constant and equal to |u;| cosh((ii), 

(iv) the hyperbolic center of mass moves along the geodesic that is perpendicular to the geodesic containing mi, 711 - 2 . 
Its velocity vector v has constant hyperbolic norm equal to |a;|. 


Proof. It is straightforward to see that (i) holds since the action is by isometries. To prove the other statements we 
work on the model. By Proposition 5.4 it is sufficient to show that the statements hold for the relative equilibrium 
in canonical configuration. Using (|4.3[) it is seen that such relative equilibrium is given by 


{xi{t),yi{t)) = e*‘^(cos 6 »i,sin 6 »i), {x 2 {f),y 2 {t)) = e“^(-cos 6 » 2 , sinda), 


where w is given by ( |5.6[ ). 

At the time t, the masses lie along the geodesic it that consists of the upper half of the circle of radius centered 
at (0,0). The velocity vector of each particle points along the radial direction and is clearly perpendicular to this 
geodesic (we are using that the model is conformal). Moreover, since e^^ sin 6*1 and sin 62 have the same sign, 
the velocity vectors are seen to have the same orientation so we have proved (ii). 

Next, the vector vi = a;e*‘^(cosdi,sin0i) and its hyperbolic norm is given by 


iyiitW 
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Substituting 2.3 into (5.6) yields the expression for w given in the statement of the theorem. Equation (2.3) also 
implies that = cosh(di). So we have shown that (iii) holds for vi. An analogous calculation shows that it holds 
for V 2 - 


Finally, note that the hyperbolic center of mass is located at e*‘^(0,1) = (0, e*“) at time t. It therefore traverses a 
vertical geodesic that is always perpendicular to it- The hyperbolic norm of its velocity vector is readily computed to 
be |a;|. 


□ 


[51 


A schematic description of the initial conditions that give rise to a hyperbolic relative equilibrium is given in Figure 



Figure 5. Schematic representation of the initial conditions for a hyperbolic relative equilibrium. The 
vectors vi,V 2 are perpendicular to the geodesic passing through mi, m 2 - Their respective hyperbolic 


norms are |a;| cosh(di) where w satisfies (7.1). During the motion, the hyperbolic center of mass moves 


along a geodesic that is perpendicular to the one connecting mi and m 2 and its velocity vector v has 
constant hyperbolic norm |w|. 


The above theorem has a nice geometrical interpretation. Consider two particles in whose velocities are per¬ 
pendicular to the geodesic that contains them and point in the same direction. There are two competing effects that 
act on the particles. On the one hand, due to their inertia, they tend to follow geodesic curves. These curves spread 
apart since we are in hyperbolic space. On the other hand, there is the gravitational effect that pulls the particles 
together. The theorem tells us that there exist unique, constant values of the particles’ speeds that exactly balance the 
above effects, and the distance between the particles remains constant throughout the motion. This type of motion is 
impossible in euclidean space since in this space parallel lines do not spread apart. 

Our discussion in Section [O] shows that all of these relative equilibria are unstable. 

The intrinsic description of elliptic relative equilibria is given by the following. 

Theorem 7.2. Consider two masses mi and m 2 in the abstract hyperbolic two-dimensional space An elliptic 
relative equilibrium corresponds to the solution characterized in the following way. 

Let di, i = 1,2 denote the distance between the mass mi and the hyperbolic center of mass of mi and m 2 , d := di-\-d 2 
the distance between the masses, and let uP' satisfy: 

2kmi 2km2 


P = 


sinh^((i) sinh( 2 d 2 ) sinh^((i) sinh(2di) 


(7.1) 


Then, throughout the motion. 
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(i) di is constant, 

(ii) the velocity vectors vi and V 2 are perpendicular to the unique geodesic that contains mi and m 2 and have 
opposite orientations, 

(iii) the hyperbolic norm of Vi is constant and equal to |w| sinh(di), 

(iv) the hyperbolic center of mass is fixed during the motion, 

(v) the motion is periodic with period . 


Proof. We proceed in analogy with the proof of Theorem 7.1 and work with a canonical conhguration in the H^-model. 

(— cos 02 cos{tuj), sin 62 ) 


According to (4.4), such relative equilibrium is given by 

(cos 9i cos{tu}), sin 61 ) 


{xiit),yi{t)) = 


1 + cos 01 sin(tu;) 


ix2it),y2it)) = 


1 — cos 02 sin(tw) 


where uj is given by ( |5.6| ). Recall that substitution of (2.3) into ( |5.6[ ) yields the expression for w given in the statement 
of the theorem. The above formulae prove (v). 


As before, statement (i) follows since the action is by isometries. Using ( |4.4[ ) with x = 0,y = 1, shows that the 
hyperbolic center of mass is fixed during the motion so we have shown (iv). 

Next, the geodesic that contains mi and m2 at time t is either the semi-circle 


it = 

or the vertical line 


{ix,y) 


: (x -I- tan(ta;))^ + y^ = sec^{tu 


it = {{x,y) : X = 0 } if 


t = 


if 

{k -t- i)7r 


{k + i)7r 


t ^ ^, fc G Z, 


to 


, fc G Z. 


In either case, it is conveniently parametrized by hyperbolic arc-length as 


tanh(s) cos(tw) 


~ l+l^h^^^sinM’ ~ 1 


sech(s) 


s G 


Using (2.3) one checks that such parametrization satisfies 
i.Xg{di),yg{di)) = {xi{t),yi{t)), 

One computes 


Xg{s) = 


sech^(s) cos(tuj) 


-I- tanh(s) sin(tw) ’ 

[xg{-d2),yg{-d2)) = {X2{t),y2{t)). 
tanh(s)sech(s) -I- sech(s) sin(ta;) 


(7.2) 


y'g{s) = - 


® ( 1 -I- tanh(s) sin(ta;))^ ’ ® (1 -I- tanh(s) sin(ta;))^ ’ 

where ' = Using ( |2.3| ) it is seen that the velocity vectors vt = {xi{t), yiit)), i = 1,2 can be written as 

cos 02 


vi = {y'g{di), -x'g{di)) , 


V2 = -OJ- 


{y'g{-d2),-x'g{-d2)) . 


sin 01 ^ ® 3 J sin 02 

Since is conformal, this shows that both vectors are perpendicular to it and are oppositely oriented (note that 
x'g{di) and x'g{—d 2 ) have the same sign). This proves (ii). Moreover, from the above expression we have that the 
hyperbolic norm of vi is 

cos 01 


w 


sin 01 


Iy'gidi)^ + x'gidi)"^ I I cos 01 

= = |cl;| sinh(di) 


Vgidiy 


sm til 


where we have used ( |2.3[ ) and that the parametrization (7.2) is by hyperbolic arc-length. Analogously, the hyperbolic 
norm of V 2 is |a;| sinh(d 2 ) and we have shown (iii). 


□ 


E 


A schematic description of the initial conditions that give rise to an elliptic relative equilibrium is given in Figure 

The stability results for elliptic relative equilibria can also be formulated in an intrinsic way. For this matter recall 
that the stability condition can be expressed by equation (6.10). Using (2.3) we conclude that the elliptic relative 
equilibrium described in Theorem 7.2 is G^-stable if 


mi 

m 2 


< 


y^3tanh^(di) -|- 1 

4sinh3(di) ’ 
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Figure 6 . Schematic representation of the initial conditions for an elliptic relative equilibrium. The 
vectors vi,V 2 are perpendicular to the geodesic passing through mi, m 2 - Their respective hyperbolic 
norms are |a;| sinh(fii) where uj satisfies (7.1). During the motion, the hyperbolic center of mass is at 
rest. 


and it is unstable if the reverse inequality holds. 

Figure [TI illustrates the values of di for which the solution is stable for a given ratio In essence, the stability is 
lost when the masses are far away. 



Figure 7. 


Stability properties of elliptic relative equilibria as a function of di and the mass ratio —. 


Remark 7.3. It is remarkable that the hyperbolic center of mass in both types of relative equilibria has a behavior 
that is similar to its euclidean counterpart. Namely, it follows a geodesic at constant speed or it stays fixed during 
the motion. As numerical calculations show, this kind of behavior should not be expected for generic solutions of the 
problem. This is usually referred to as the “absence of the integral of center of mass” (see [TT]b 

o 


Appendix. 

In this Appendix, we will prove some general results on G^-stable relative equilibria for symmetric Hamiltonian 
systems. Although these results seem to be quite standard, we include them for the sake of completeness. 
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Let (P, ri) be a symplectic manifold with symplectic 2-form G a Lie group which acts symplectically on P, 
J : P —>■ g* an equivariant momentum map and H : P —>■ K a G-invariant Hamiltonian function. 

Suppose that z : K —>■ P is a relative equilibrium, 

z{t) = exp(^t) • z(0), with ^ e g, 

and that J(z(0)) = fi. Using Noether Theorem, the equivariance of the momentum map J and proceeding as at the 
beginning of Section]^ we deduce that z(IR) C G^ • z(0), that is, z(]R) is contained in the G^-orbit by z(0). 

Next, assume that z is a G^-stable relative equilibrium. This means that for every G^-invariant open subset U of 
P, with z(]R) C P, there exists an open subset W oi P such that z(M) C lU C P and for every solution z' of the 
Hamilton equations with initial condition in W, we have that z'{t) G P, for all t (see [5T|L 
Now, we will consider two special cases: 

1. The isotropy group G^ is compact. In this case, we will prove that solutions sufficiently close to z are 
bounded. In fact, we will see that the following result holds. 

Proposition A.l. If z :M. ^ P is a Gp,-stable relative equilibrium and Gp is compact then there exists an open subset 
W ofT*Q, such that z(]R) C W, and for every solution z' of the Hamilton equations with initial condition in W we 
have that z' is bounded. 

Proof. Let Pq be an open subset of P, with z(0) S Pq, and K a compact subset such that Pq C K. Consider the 
subset U oi P given by 

U = Gp- Po. 

It is clear that P G Gp ■ K and Gp ■ K is compact. Thus, P is bounded. Moreover, P is G^-invariant and z(IR) 
Therefore, there exists an open subset W in P, such that for every solution z' with initial condition in IT, we 
that z'{t) G P, for all t. Now, using that P is bounded, we conclude that z' also is a bounded curve. 

2. The Hamiltonian system is the 2-body problem in an arbitrary Riemannian manifold. In this 
P is the cotangent bundle of a manifold Q endowed with the canonical symplectic structure, where 

Q = (M X M) \ A, 

with M an arbitrary Riemannian manifold and A the diagonal in M x M . 

Moreover, we will suppose that G is a Lie subgroup of the isometry Lie group of M and the action of G on T*Q is 
the cotangent lift of the diagonal action of G on Q. 

In addition, the Hamiltonian function H : T*(3 —>■ K on T*Q is given by 

H{pi,P 2) = ^\\pi\\^ + ^\\p2f + V{d{XuX2)), 

Zttii 2 tti 2 

with mi, m 2 > 0, pi G Tf^M, p 2 G Tf^M, {xi,X 2 ) G Q and V : K+ —>• M a smooth function. As in the paper, || - || 
denotes the norm in the fibers of T*Q induced by the Riemannian metric. On the other hand, d : Q —> K is the 
restriction to Q of the Riemannian distance in M x M. 

Now, denote by r the map defined hy r = do ttq, with ttq : T*Q ^ Q the canonical projection. Then, since z(M) 
is contained in the G^-orbit of z(0) and G is an isometry group, we have that 

(r o z){t) = (r o z)(0 ) := Tq > 0 , for every t S M. 

Furthermore, we may prove that the Riemannian distance between two particles along a solution sufficiently close to 
the G^-stable relative equilibrium z is almost equal to rg. 

Proposition A.2. For every e > 0, there exists an open subset W in T*Q such that z(K) C W and for every solution 
z' of the Hamiltonian system {T*Q,H) with initial condition in W we have that 

|r(z'(t)) - rol < e, Vt. 


open 


C P. 
have 
□ 

case. 


Proof. Consider the open subset of T*Q 


U = r ^{ro- e,ro + e). 
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Then, U is G-invariant and z(M) C {/. As a consequence of Gy^-stability of the relative equilibrium z, there exists an 
open subset W CT*Q, with z(IR) C W, such that 

z'{t) G U, Vt, 

for any solution z' with initial condition in W. Therefore, 

|r(z'(t)) - rol < e, Vt. 

□ 

It is easy to prove that a natural extension of the previous result also holds for the more general case of the n-body 
problem in M, with n > 2. 
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